Non-empirical pairing energy functional in nuclear matter and finite nuclei 
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We study 1 So pairing gaps in neutron and nuclear matter as well as T = 1 pairing in finite nuclei 
on the basis of microscopic two-nucleon interactions. Special attention is paid to the consistency 
of the pairing interaction and normal self-energy contributions. We find that pairing gaps obtained 
from low-momentum interactions depend only weakly on approximation schemes for the normal 
self-energy, required in present energy-density functional calculations, while pairing gaps from hard 
potentials are very sensitive to the effective-mass approximation scheme. 
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I. INTRODUCTION 



Medium-mass and heavy nuclei can be studied system- 
atically through nuclear energy-density-functional (EDF) 
calculations [1]. Within a single-reference implementa- 
tion, the minimization of the energy-density functional 
leads to solving Hartree-Fock-Bogoliubov (HFB) equa- 
tions 0. However, nuclear energy functionals accounting 
for correlated single-particle motions and superfluidity 
employed so far are of (semi-)empirical character [![. It is 
a central goal to construct non-empirical energy-density 
functionals connected to two- and many-nucleon inter- 
actions in free space Q in view of the challenges posed 
by exotic nuclei with an unusually large ratio of neu- 
trons over protons. The development of low- momentum 
interactions based on renormalization group (RG) meth- 
ods [1, []| opens up such a possibility, as they enable 
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technically simpler many-body approaches 

Such a long-term project to connect the nuclear EDF 
to underlying nuclear interactions has recently been ini- 
tiated, first focusing on the part of the energy functional 
that drives pairing properties of nuclei [H, E3, EH- A 
difficulty is that a quantitative description of superflu- 
idity in nuclear systems is a delicate task that a pri- 
ori requires the treatment of complicated many-body 
processes. In fact, an on-going discussion concerns the 
impact of medium polarization effects, beyond the di- 
rect term of the nucleon-nucleon (NN) interaction, on 
pairing properties of finite nuclei. In Refs. [12, EH, 
two thirds of the observed neutron pairing gap was ac- 



counted for in 120 Sn by using the Argonne V14 NN po- 
tential (T3 | as pairing interaction and combining this 
with the semi-empirical Skyrme functional SLy4 in the 
particle-hole channel. Adding induced interactions and 
self-energy effects due to the exchange of collective fluctu- 
ations between nucleons moving in time-reversed states, 
the missing one third was recovered [H, EH, EH Ell- 
in Refs. (IP . EH, however, neutron and proton pair- 
ing gaps were found to be consistent with experimental 
data over a large range of semi-magic nuclei when low- 
momentum NN interactions Mow k were used as pairing 
interaction, combined with the same Skyrme functional 
in the particle-hole channel. In this case, one therefore 
expects that neglected many-body forces and polariza- 
tion effects result in a small net contribution to pairing 
gaps in (known) finite nuclei. 1 

Before addressing the contribution of many-body 
forces and collective fluctuations to pairing gaps, the aim 
of the present work is to understand the qualitative and 
quantitative mismatch between the two sets of results 
published in Refs. [13, EH and Refs. [H,EI|], respectively, 
which both employ a free-space NN interaction as pairing 
force and the semi-empirical Skyrme functional SLy4 in 
the particle-hole channel. The only difference between 
the two calculations resides in the intrinsic resolution 
scale (in the RG sense [1, [Bj]) of the Argonne V14 and 
Mow k NN interactions, which both reproduce the relevant 
low-energy scattering phase shifts [17] . Starting from 
this observation, one is led to focus on the coupling of 
the pairing interaction to the normal self-energy, in par- 
ticular when the latter is approximated by a momentum- 
independent effective mass as is the case for present EDF 
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calculations. Settling this issue requires fully microscopic 
calculations, where both the normal and anomalous self- 
energies are computed consistently from the NN inter- 
action. To date, this is only possible for infinite nuclear 
matter (INM) (see, for example, Ref. [l8j). Of course, the 
results obtained in INM cannot be extrapolated straight- 
forwardly to finite nuclei. Nevertheless, calculations of 
pairing gaps in INM provide a base-line for the schemes 
used in finite nuclei in Refs. 

nana 

, [HI, EH and allow us 
to probe the sensitivity of pairing gaps to approximations 
of the normal self-energy. 

We stress that our goal is not to perform the most in- 
volved calculations of pairing gaps in INM, for example, 
including induced interaction and associated self-energy 
effects [ill HO, [U [22[. Rather, we work at lowest order 
in the many-body expansion, with special attention (i) 
to the consistency between the normal and anomalous 
self-energies when using either a hard or low-momentum 
NN interaction and (ii) to the effects of neglecting the 
momentum dependence of the effective mass and quasi- 
particle strength when solving the gap equation. 

This paper is organized as follows. In Section [Til the 
many-body frameworks used to expand the normal self- 
energy and pairing gap are set up for NN interactions 
characterized by either a low- or high-momentum resolu- 
tion scale (in the RG sense [3, Q). Procedures to aver- 
age the momentum dependence of the normal self-energy 
are discussed. Section IIIII presents results for the effec- 
tive masses and pairing gaps in pure neutron matter and 
symmetric nuclear matter. The impact of the momen- 
tum averaging of the normal self-energy on pairing gaps 
is analyzed. In Section IIV1 we discuss the consequences 
of our findings in INM on the computation of neutron 
and proton gaps in semi-magic nuclei. We conclude and 
give an outlook in Section [V] 



II. MANY-BODY FRAMEWORK 
A. Hamiltonian 

The basic ingredient to microscopic calculations is the 
Hamiltonian that incorporates two- and many-nucleon 
interactions constrained by scattering experiments and 
few-body properties. In the present work, we neglect 
many-nucleon interactions, although it is important to 
characterize their impact on pairing properties [23| . 
Many-body forces may change the value of the gap, es- 
pecially toward higher density, but are not expected to 
alter the conclusions of this paper. 

The setup of a meaningful expansion scheme for nu- 
clear many-body calculations depends on the choice of 
NN interaction [6j, [7], |8[ . We consider two schemes that 
are currently used in low-energy nuclear structure and 
reaction calculations. First, we work in a scheme which 
attempts to model the short-range parts of nuclear forces 
explicitly and is thus characterized by a large intrinsic 
resolution scale Ah ar d- This will be referred to as a "hard" 



NN interaction. Second, we consider low-momentum 
interactions with lower intrinsic resolution scale A so f t , 
which we refer to as "soft" NN interactions. 

For both cases, we generate the interaction matrix ele- 
ments starting from the Argonne via potential [24j by 
solving the symmetrized RG equation [f| [TtJ with an 
exponential regulator /(fc,A) = exp [—(&;/ A) 2riexp ] with 
n cxp — 7. The hard interaction is obtained by evolving 
to Ahard = 6.0 fm -1 . Using such a cutoff-scale instead of 
the initial Argonne vi% potential allows one to reduce the 
numerical complexity of the calculations while maintain- 
ing the features of a hard potential. The soft interaction 
Mow k is obtained by further evolving the RG equation to 
a typical scale A so f t = 1.8 fm -1 [3,0- 



B. Expansion scheme 

Starting from hard interactions, the short-range parts 
must be summed before the interaction can be used in 
many-body calculations. The traditional approach is 
based on the Brueckner G matrix [25| and on the reorga- 
nization of binding- or self-energy expansion schemes in 
terms of the number of hole lines entering the retained 
diagrams '26, 27[. In addition, there are calculations for 
nuclear matter and finite nuclei based on a self-consistent 
in-medium T matrix, which sums particle-particle as well 
as hole- hole ladders (see, for example, Ref. [28|). How- 
ever, both expansions in terms of the G or T matrix are 
nonperturbative for hard interactions, because they do 
not decouple low and high momenta [6|, |7( . 

Starting from low-momentum interactions [3, Q the 
high-momentum modes are decoupled. This offers the 
possibility to use a perturbative expansion, as was shown 
explicitly for the particle-particle-channel contributions 
to the energy Q. It remains to be checked that the ex- 
pansion is perturbative in the particle-hole channels. 

It is thus crucial to realize that the many-body expan- 
sion scheme differs depending on the intrinsic resolution 
scale characterizing the Hamiltonian. A trivial, but es- 
sential, implication is that only the complete resumma- 
tion expressed in terms of the full Hamiltonian provides 
results that are independent of the expansion scheme, 
whereas results computed at a given order in the relevant 
expansion may differ depending on the scheme used. This 
underlines the necessity to specify the scheme employed 
and to perform consistent calculations, such as comput- 
ing normal and superfluid self-energies at the same order 
in the relevant expansion scheme. In addition, due to 
the necessity to rearrange the expansion scheme depend- 
ing on the resolution scale A, results obtained through 
truncated calculations cannot be expected to be cutoff 
independent. 

Table fl] compares the set of diagrams taken into ac- 
count to compute the normal and anomalous self-energies 
to first and second order in the two expansion schemes 
considered. In both cases, the first-order anomalous di- 
agram is characterized by the use of the direct NN in- 
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Hard NN interaction 
(large A) 



s£L = Ho 
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TABLE I: Expansion schemes for soft (left) and hard (right) 
interactions to first (up) and second (down) order. The 
dashed line denotes the free-space NN interaction and the 
wiggly line the G matrix. Diagrams with more than one 
anomalous propagator are not shown. 



teraction as the (particle-particle-irreducible) kernel en- 
tering the gap equation. The first-order normal self- 
energy is provided by the Hartree-Fock (HF) diagram 
when starting from soft interactions, whereas it leads to 
the Brueckner-Hartree-Fock (BHF) approximation when 
employing hard potentials. The second-order diagrams 
correspond to a non-collective treatment of screening and 
vertex corrections as is done in Ref . [22| . 

The goal of this paper is to study 1 So pairing gaps in 
nuclear matter and finite nuclei at lowest order in the two 



expansion schemes, with special attention to the consis- 
tency of the pairing interaction and normal self-energy 
contributions. In particular, we focus on the effect of the 
normal self-energy and of effective-mass approximation 
schemes on the pairing gaps. 



C. Normal self-energy 

At lowest order, the normal self-energy takes the form 
sST(Pi)= $> T '(p 2 )(k|^'|k>, (1) 

P2.1"',o-' 

s££(Pi,")= E" T W(k|G-'(P,c + e ^ 2 ))|k) , 



P2,T',<7' 



(2) 



with the notation p = |p| and where the occupation func- 
tion n T (p) = 8 (hp — |p|) is taken as a step function. 
The relative and center-of-mass momenta are defined by 
k = (pi — p 2 )/2 and P = pi + p 2 , respectively. The in- 
dex r characterizes the single-particle isospin projection. 
Using normal-state Fermi-Dirac distributions functions 
rather than BCS-like occupation numbers is a satisfac- 
tory approximation in INM around saturation density 
where the pairing gap A T is small compared to the Fermi 
energy ef = {k T F ) 2 /{2m) j^. 

The Brueckner G matrix is calculated through a 
partial- wave expansion (in units of h = c = 1) 



J 



{k'\ G$ SJ {P,w) \k) = (k'\ VSTsj \k) + - <l 2 dq V (fc'| V%, \q) 



i 



r 



where the angular-averaged Pauli-blocking operator and 
two-particle-state energies are defined as 



{Q^'{P,k)) = 



1 



G?COS ( 



'Pk 



l-n T (pi) l-n T (p 2 ) 



(e TT '(P,k)) = - dcosflpk [e T (pi) +e T \pi) 



The on-shell single-particle energy e T (p) entering Eqs. (|2|) 
and ^ is obtained through 



to pure neutron matter (PNM) and symmetric nuclear 
matter (SNM) in the present work. Because we focus 
on neutron-neutron pairing, only neutron self-energies, 
r = n, are eventually needed. Thus, the short-hand nota- 
tions e(p) = e"(p),E HF (p) = and S BHF (p,w) = 

S&(|)| ijj) are used, along with corresponding notations 
for effective masses, quasiparticle strength, pairing gaps 
and Fermi momenta introduced below. 



e T (p) 



2m 



Re £«>,£»). 



(4) 



Equations ©, © and ^ are solved self-consistcntly 
when using hard interactions. For soft interactions, the 
system reduces to the direct evaluation of Eqs. |T|) and (j4|) 
thanks to the energy independence of S^^". 

We note that these equations are valid for asym- 
metric nuclear matter, although we restrict ourselves 



D. Effective- mass approximation 

1. Definitions 

A focus of the present work is to study effective- 
mass approximation schemes on pairing gaps. The 
momentum-dependent effective mass m* (p, k F ) is de- 
fined by 



4 



de T (p) _ p 
dp 



m r(P: fc F) ' 



(•5) 



surface. The free single-particle spectrum £g(p) = (p 2 ~ 
fcp 2 )/(2m) and a typical width of A = 2.0 MeV are used 
for simplicity in Eq. (|10p . 



where fcj denotes the dependence on both Fermi mo- 
menta. This total effective mass can be separated into 
the product of the fc-mass and the e-mass defined by (30l | 



m r,k(P> fc F ) 

m 

™*,efe fc F ) 



l+— ^eT, T A (p,uj) 



= 1 - 



p dp 



, (6) 
(7) 



'(p) 



The fc-mass relates to the spatial non-locality of the 
normal self-energy, whereas the e-mass characterizes the 
dynamical correlations associated with the energy depen- 
dence. The e-mass can also be expressed in terms of the 
quasiparticle strength, or Z-factor, 



Z T (p, fcp 



(8) 



which quantifies the quasiparticle part of the one-body 
Green's function and occurs in the pole approximation 
discussed in the following. 



2. Momentum-independent approximation 

At this point, the introduction of the effective mass 
is essentially a matter of definition. The real purpose 
is usually to neglect its momentum dependence in or- 
der to recover a (density-dependent) quadratic disper- 
sion relation for e T (p). In the present work, it is mo- 
tivated by the need to make the connection with the 
Skyrme functional whose self-energy can at best be re- 
lated to a momentum-independent approximation of the 
microscopically-obtained effective mass. 

Hence we have to reduce the momentum dependence of 
m*(p, kp ) and Z T (p, fcp ) to a dependence on fcp. Obvi- 
ously there are different ways to do so. In order to probe 
the sensitivity of observables to the particular scheme 
used, we consider two approximations. The first, which 
we denote as the point-evaluation ( "pe" ) approximation, 
is standard and consists of taking the value at the Fermi 
momentum, with X T = m* or Z T , 



Xp e (kp ) — X T (p — fcp, fcp ) . 



(9) 



The second method proceeds through an averaging 
( "aw" ) of the momentum dependence over the Fermi sur- 
face, 

a " (F) - J f(q t A) ^dq iirvr ' (10) 



where u T q v T q = A/(2y / (Q(q)) 2 + A 2 ) denotes the BCS 
pair occupation function, which is peaked at the Fermi 



E. Anomalous self-energy 

After discussing the computation of the normal part of 
the self-energy together with its effective-mass approxi- 
mation, we now turn to the anomalous self-energy that 
leads to pairing gaps. We are interested in neutron- 
neutron pairing at sub-saturation densities fcF ~ 0.8 — 
1.4 far 1 in PNM and SNM. Neutron superfluidity in this 
density range in INM is in the % channel [13, EH- In 
this work the pairing kernel is restricted to the direct 
NN interaction and therefore we keep only the 1 So par- 
tial wave as pairing interaction. 

The gap equation is solved within the pole approxi- 
mation, which provides a good approximation to the so- 
lution of the full off-shell gap equation when the mo- 
mentum dependence of the effective mass and of the Z- 
factor is taken into account [32|. Furthermore, and as 
already mentioned, the normal self-energy and quasipar- 
ticle strength are computed in the normal state, which is 
valid for the density range considered where A/ep <C 1. 
In summary, the neutron anomalous self-energy and thus 
the gap A(p) is the solution of [33| 



AW — 1/ 



dq q' 



(11) 



with 



2m M 2 



ReE A (p, £ (p))+Re£ A (p, 2/x-e(p)) 



and the chemical potential is defined by /i = e(kp), so 
that it includes the normal self-energy shift with respect 
to the free Fermi energy. The chemical potential can also 
be calculated self-consistently to account for the effect of 
pairing correlations, but we have checked that this has 
very little effect in the density range of interest. 

It is important to realize that for Z(q) ^ 1 the physical 
gap A (p = fcp) in the excitation spectrum of the system 
is given by [34| 



A(p) = Zip) A(p) . 



(12) 



Linearizing additionally Sa (p, cj) in energy around \x 
leads to the BCS-type equation 

fi M 1 f, 2 Z *(P) (p|^k)^(a)A,( g ) 

Ai(p) = / dqq = , (13) 

nJ Vtf(?)+a?(g) 

where = £i{p) - p. 

The index i in Eq. (JT3J) labels different cases considered 
in this paper regarding the choice of single-particle en- 
ergy and quasiparticle strength. We define three classes 
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(i) Asoft 


0(p)=£ HF (p)-M 

&/s(p) = (P 2 ~ k 2 F )/(2m;?f av (k F )) 


Zi(p) = 1 
Z2/3 (p) = 1 


(ii) A;f ard 


£i(p) = e BHF (p)-/i 

6/3(p) = (p 2 -fc F )/(2m;^ F „(fc P )) 


^( P ) = Z BHF ( P ) 
Z 2/3 (p) = Z™ F av (k F ) 


(iii) ALrd 


t2 /3 (p) = (p 2 -kl)/(2ml B «J av (k F )) 


Zi( P ) = 1 

^2/3 (p) = 1 



TABLE II: Three classes of calculations: (i) first order in a 
soft NN interaction, (ii) first order in a hard NN interaction, 
(iii) first order in a hard NN interaction but neglecting all 
effects related to the energy dependence of the normal self- 
energy. Each class contains a calculation of reference retaining 
the full momentum dependence of the normal self-energy (i = 
1) and the pe/av schemes (i = 2, 3). 



depending on whether one starts from a soft or a hard 
NN interaction, see Table UT1 For each class, one goal is to 
compare the gaps obtained from strictly solving Eq. (fT3|) 
to those obtained using further approximations, for ex- 
ample neglecting the momentum dependence of the ef- 
fective mass and of the Z-factor. The third class (A 1 ard ) 
defined in Table |TT] is not consistent in the sense that we 
only keep the fc-mass of the total BHF effective mass and 
neglect all e-mass and Z-factor effects in the gap equa- 
tion 2 . This can be considered as an intermediate case 
between classes A* ft and K? rA . 



III. INFINITE NUCLEAR MATTER RESULTS 
A. Effective mass 

1. Momentum dependence 

Figure [1] shows the momentum dependence of the ef- 
fective masses calculated at lowest order for soft (upper 
panel) and hard (lower panel) interactions, at a represen- 
tative density of fc F = 1.2 fm" 1 in PNM and SNM. 

The HF calculation using the soft interaction only 
generates a fc-mass that is smaller than the bare mass 
and displays a smooth momentum dependence. The 
HF fc-mass is smaller in SNM than in PNM due to 



2 Once the fc-mass is extracted from the total effective mass, the 
corresponding self-energy is recovered through the integral given 
in Table |II] 
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/ \ e-mass, SNM - 
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t e-mass, PNM - 
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FIG. 1: (Color online) Momentum dependence of the HF ef- 
fective fc-mass for A so f t = 1.8 fm -1 (upper panel) and of the 
BHF effective masses for Ahard = 6.0 fm -1 (lower panel) in 
PNM and SNM at k F = 1.2 fm" 1 . Whereas the HF approx- 
imation only provides a fc-mass, the hole-line expansion al- 
ready generates an e-mass at leading order. 

the stronger in-medium effects, notably brought by the 
proton-neutron (tensor) interaction. 

In the BHF calculation already at leading order a fc- 
mass and an e-mass are generated. The e-mass is given 
by the energy dependence of Re s[ i ^. d (p, w) and displays 
a typical enhancement around the Fermi momentum as- 
sociated with the increased probability of virtually oc- 
cupy two-particle-one-hole configurations. This effect 
comes at second order in the perturbative expansion for 
soft interactions (see Table U}. The BHF fc-mass is sim- 
ilar to, but slightly smaller than the HF fc-mass shown 
on the upper panel. The total BHF effective mass is the 
product of the fc-mass and the e-mass and is thus larger 
than the fc-mass for all densities, in addition to carrying 
the typical enhancement of the e-mass around the Fermi 
momentum. We have also checked that BHF masses and 
Z-factors do not change significantly by increasing the 
cutoff beyond A = 6.0 fm" 1 in SNM and PNM. 

2. Averaged momentum dependence 

Figure [5] shows the momentum-independent effective 
masses Trit^Qtp) and (fcp) obtained by applying 
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FIG. 2: (Color online) Momentum-independent effective 
masses TOpJ 1 F (fcp) and m^ F (fcp) obtained from the soft in- 
teraction in PNM and SNM. For comparison we show the 
effective masses mg Ly4 of the Skyrme SLy4 parameterization. 

the two averaging schemes introduced in Section Ul D 2l for 
the soft interaction. A key result obtained for PNM and 
SNM is that the two schemes lead to essentially identical 
results. This is due to the mild momentum dependence 
of the HF fc-mass obtained from the soft interaction and 
indicates that reducing such a momentum dependence 
may be a tractable approximation when solving the gap 
equation. Figure [2] also compares the microscopically- 
calculated 77i* HF (fcp) to the density-dependent effective 
mass ^sly^^f) 01 ^ ne SLy4 Skyrme functional param- 
eterization. Although the functional dependence is not 
fully captured by TOsLy4(^F), the microscopic result in 
SNM at saturation density fcp = 1.35 fin - is well re- 
produced. In PNM, however, the neutron effective mass 
m SLy4(^f) underestimates the microscopic predictions 
significantly. This reflects the known deficiency in the 
isovector effective mass of SLy4 35] . 

In Fig. [21 we present the momentum- independent total 
effective mass, fc-mass, and Z-factor obtained in SNM ap- 
plying the two averaging schemes in the case of the hard 
interaction. Due to the more pronounced momentum de- 
pendence of the effective masses and the larger averaging 
region set by the regulator f(q, A = 6.0 fm _1 ) in Eq. (flU)) . 
the point-evaluated and averaged values differ substan- 
tially. This difference is also much larger in SNM than in 
PNM as self-energy effects are larger in SNM. Due to a 
compensation effect between the fc-mass and the e-mass, 
the total effective mass happens to be relatively insensi- 
tive to the averaging scheme used. Nevertheless, as the 
gap equation depends on both the total effective mass 
and Z-factor (see Eq. (|T3"]) ). the present results indicate 
that it may be unreliable to average the momentum de- 



FIG. 3: (Color online) Momentum-independent effective 
masses and Z-factors obtained from the hard interaction in 
the BHF scheme in SNM. m* k mF (fa?) and Z BHF (k F ) depend 
strongly on the averaging scheme, whereas the total mass is 
less sensitive to this. 

pendence of m* BHF (p, kp) and Z(j>, kp) in this case. The 
same conclusion can be anticipated when solving the gap 
equation within the frame of the third class of Table [Til 
that is when the e-mass and Z-factor are neglected. 

B. Pairing gaps 

1. Soft interaction 

Figure Q] shows our results for the pairing gaps at the 
Fermi surface A(p — kp) as a function of the Fermi mo- 
mentum kp in PNM and SNM for the Ag oft class de- 
fined in Table [TTJ The band represents the pairing gaps 
obtained using the SLy4 effective mass and varying the 
cutoff in the NN interaction over a wide cutoff range. 
This band corresponds to the variation of the pairing 
gaps obtained in finite nuclei (36| by varying the resolu- 
tion scale of the pairing kernel. The upper limit of the 
band corresponds to the low cutoff A = 1.8 fm -1 . This is 
the result to be compared to the microscopic calculations 
discussed in this subsection. The lower limit of the band 
corresponds to the hard cutoff A = 15.0 fm _1 and will be 
relevant to the next subsection. 

As shown in the upper panel of Fig. [H the pairing gaps 
obtained from the Ag oft class are essentially indistinguish- 
able. We find 

A[e HF (fc)] « A[m;« F (M] « A>:« F (M] , (14) 

to an excellent approximation in PNM. As expected from 
Fig. [21 the momentum averaging of the normal self- 
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0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 
k ¥ [fm- 1 ] 

FIG. 4: (Color online) Neutron 1 So pairing gaps in PNM and 
SNM obtained using the soft interaction with A so f t = 1.8 fm -1 
as the pairing interaction and the HF approximation for the 
normal self-energy. Results are shown for the three cases of 
class A soft (see Table ITT)) . 



energy e HF (p) — > (p 2 — fcp)/(2m* HF (/cp)) is well justified 
in PNM and has essentially no impact on the computed 
gaps. Due to the wrong isovector dependence of mg Ly4 
compared to m* HF (fcp) (see Fig. the microscopic gaps 
are larger than the upper limit of the band. 

In SNM the pairing gaps are also insensitive to the 
effective-mass approximation scheme, in particular over 
the range kp « 1.0 — 1.4 fm - . In the density region 
kp w 1.2 — 1.4 fm - , TigL y 4 reproduces well the effective 
mass obtained from the soft NN interaction (see Fig. [2]) 
and therefore the microscopic gaps are close to the upper 
limit of the band. At lower densities kp < 1.2 fm - , 
TO SLy4 i s larger than the calculated effective mass and 
therefore we find pairing gaps that are smaller than the 
upper limit of the band. 

In addition, we practically find cutoff independence 
of the pairing gaps in SNM and PNM for soft cutoffs 
A w 1.8 — 3.0 fm - . This approach therefore provides a 
tractable lowest-order starting point with respect to ap- 
propriate 3 variations of the renormalization scale. These 



_ i i i | i i i | i i i i i | i i r^T i i | i i i | i i i _ 




0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 
k F [fm" 1 ] 

FIG. 5: (Color online) Neutron 1 So pairing gaps in SNM ob- 
tained using the hard interaction with Ah ar d = 6.0 fm -1 as 
the pairing interaction and the BHF approximation for the 
normal self-energy. For the results in the upper panel, only 
the fc-mass effects are taken into account according to class 
Ah ar d of Table [TTJ The lower panel shows the results of the 
class A{f ard including e-mass and Z-factor effects. 

results also complement the cutoff independence of pair- 
ing gaps obtained in INM using a free single-particle spec- 
trum |l7l ] and in finite nuclei [3(| over the same cutoff 
range. 



2. Hard interaction 

Figure [5] shows the pairing gaps A(p = kp) computed 
in SNM according to the two classes A^ ard and Ajf ard de- 
fined in Table HT1 The upper panel shows the results tak- 
ing only the fc-mass contributions into account (Aj ard ), 
whereas the lower panel includes also the e-mass and Z- 
factor effects (A^ ard ) generated at first order in the hole- 
line expansion. 

Comparing the two different classes, we find the typi- 
cal systematic reduction of the gaps due to the decreased 
spectral strength of the quasiparticle propagator that 



3 One must consider variations of A such that the perturbative 



expansion remains valid. Weinberg eigenvalues demonstrate that 
this is the case for A < 3.0 fm -1 |37ll . 
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wins over the increased density of states characterizing 
the e-mass effects. At this point, however, we are not 
primarily interested in the differences of the pairing gaps 
between the different classes of Table [TTl but rather in the 
deviation of the pairing gaps within a given class as we 
approximate the momentum dependence of the normal 
self-energy and of the Z-factor. 

It is clear from Fig. [5] that the choice of the method 
used to average the momentum dependence of the ef- 
fective mass and of the Z-factor has a strong impact 
on pairing gaps, irrespective if e-mass effects are taken 
into account or not. This result could have been ex- 
pected from Fig. [3] In contrast to the soft interac- 
tion case (lower panel of Fig. |4] for SNM), replacing 
- BHF (p) -> (p 2 - fc|)/(2m* BHF (/c F )) is unreliable. Since 

,*BHF 



fcl 



£ . 

«^*? HF is systematically larger than m* BHF (see Fig. [3]), 



we find for hard potentials, in contrast to Eq. (fT4]) . 
A[r 



Kf F (fcF)] > A[ £ * BHF (fc)] > A[m; BHF (fc F )] 



(15) 



for all densities of interest. Such an inequality holds ir- 
respective if e-mass effects arc taken into account or not. 
We also observe that the pairing gaps obtained for the 
two classes are bounded from above, over the entire den- 
sity range, by the lower limit of the band obtained when 
using the single-particle spectrum generated by the SLy4 
Skyrme functional and without explicit Z-factor. 



3. Analysis 

In order to get a deeper understanding of the quali- 
tative difference between soft and hard interactions, we 
write the gap equation, Eq. (113|) 



as 



A(p = fc F 



dqY(k F ,q) . 



(16) 



Once the self-consistent gap equation is solved, the inte- 
grand function Y(k-p,p) contains information about the 
momentum scales from which the gap at the Fermi mo- 
mentum is built. The NN interaction matrix elements 
V(kF,p), the resulting function Y(Jcf,p), and the k- 
mass are shown in Fig. [S] for cutoffs A so f t = 1.8 fm -1 
and Ahard = 6.0 fm -1 at a representative density of 
fcp = 1.2 fm - . For simplicity, we only consider the class 
A^ ard for the hard interaction. This way, we avoid the ad- 
ditional subtleties connected with Z-factors and the more 
pronounced momentum dependence of the total effective 
mass for the hard potential case. 

As can be seen from the middle panel of Fig. [6j the 
gap is generated for soft interactions from momentum 
modes around the Fermi surface since off-diagonal ma- 
trix elements do not couple low and high momenta (up- 
per panel). It is therefore understandable that fixing 
the momentum-independent effective mass to its value 
at or in the vicinity of the Fermi surface (lower panel) is 
a good approximation. For large cutoffs, however, the 
major contributions to the gap at the Fermi momen- 
tum originate from high momenta [l8l |. far away from 
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FIG. 6: (Color online) Analysis of the effective-mass approx- 
imation at a representative density = 1.2 fm - . Middle 
panel: (i) for low cutoffs the gap at the Fermi momentum is 
built from low-momentum modes around the Fermi surface, 
and (ii) for large cutoffs the gap is built mainly from high- 
momentum modes, where the NN interaction matrix elements 
are maximal and repulsive (upper panel). In the hard poten- 
tial case, the momentum averaging of the effective mass would 
necessitate a non-trivial fine tuning to reproduce the pairing 
gaps calculated with the full momentum dependence of the 
fc-mass (lower panel). 



the Fermi surface. In particular for A = 6.0 fm^ 1 and 
fcp = 1.2 fm~ , the matrix elements of the pairing in- 
teraction are very small around the Fermi surface and 
essentially the entire gap strength is built from the re- 
pulsive part of the NN interaction at high momenta. 
Therefore, approximating the normal self-energy through 
a momentum-independent effective mass m*™^ defined 
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in the vicinity of the Fermi surface is unreliable in this 
case. Since the effective mass approaches the free mass m 
at high momenta (see lower panel), the gaps A[m*g HF ] 
are too small compared to the reference ones obtained 
by keeping the full momentum dependence of the nor- 
mal self-energy. In contrast, the gaps A[m*^ HF ] are too 
large compared to the reference gaps since the averag- 
ing method, Eq. (|10p. includes high- momentum modes 
with too much weight, where the effective mass is larger 
than at the Fermi surface. Hence, in order to obtain a 
more reliable momentum-independent effective mass for 
hard potentials, the behavior of the NN interaction ma- 
trix elements would have to be taken into account in the 
averaging scheme. 



IV. FINITE NUCLEI RESULTS 

Before going to higher orders, our goal is to perform a 
complete first-order calculation in superfiuid nuclei, in- 
cluding the normal and anomalous self-energies consis- 
tently. Intermediate steps towards this goal have been 
taken recently. On the one hand, calculations of non- 
superfluid nuclei have been performed starting from soft 
NN interactions [38|, l39( . This gives access to the single- 
particle field at first order in the NN interaction but 
does not account for the pairing channel. Recently, the 
work of Ref. [3£| has been extended to HFB calcula- 
tions of superfiuid nuclei, employing a soft NN inter- 
action in the particle-hole channel but without three- 
nucleon (3N) forces [40| . On the other hand, HFB cal- 
culations combining a pairing kernel based on the direct 
NN interaction with the single-particle field generated 
by the empirical SLy4 Skyrme functional have been per- 
formed [13, HH, 03, [H ED]. Eventually, one needs to 
combine the benefits of these two applications to obtain 
complete first-order calculations of superfiuid nuclei, in- 
cluding also the effects of 3N interactions on the single- 
particle field. Although the density of states around 
the Fermi energy may be dominated by NN interactions, 
3N interactions contribute especially to spin-orbit split- 
tings in nuclei. The present work is an attempt to qual- 
ify whether or not the two sets of results published in 
Refs. [13, [H| and Refs. [H, G3 provide a good approxi- 
mation to a complete first-order calculation and to under- 
stand the mismatch between them. While insights have 
been obtained in previous sections through calculations 
in INM, this present section is dedicated to assessing the 
situation in finite nuclei. 



A. Soft interaction 

Let us start with the case of the soft NN interaction. 
The findings discussed in Section [HI B H are useful to as- 
sess the validity of the results of Refs. [l(| [Hj]. One 
expects the momentum-independent effective-mass ap- 
proximation, which is reliable in INM, to be tractable 
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FIG. 7: (Color online) New Skyrme effective masses con- 
strained to reproduce the HF result obtained from the soft 
NN interaction. 
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FIG. 8: (Color online) Neutron 1 So pairing gaps in SNM and 
PNM based on the refitted Skyrme EDF (see text) with the 
soft NN interaction as the pairing kernel. 



in finite nuclei as well, because the discussion of Sec- 
tion QTTbII] can be carried over to finite nuclei. This gives 
some confidence that computing the single-particle field 
from an EDF characterized by a momentum-independent 
effective mass is a good approximation, as long as one is 
working consistently with a low resolution scale. Nev- 
ertheless, this must be checked explicitly by comparing 
the results thus obtained with the HF single-particle field 
computed in finite nuclei 4 from soft NN [38|, [39J and 3N 



4 References [38l . |39| focus on the single-particle spectrum of 
40 Ca, which displays significant differences with the correspond- 
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interactions. 



1. Refitted Skyrme effective mass 

Given the validity of the momentum-independent 
effective-mass approximation, we construct a new em- 
pirical EDF so that it reproduces well the effective mass 
obtained microscopically. As discussed, to this end the 
SLy4 value m* jm = 0.7 at nuclear saturation density 
is appropriate. However, the SLy4 effective mass has 
an incorrect isovector dependence and an unsatisfactory 
low-density behavior in SNM (see Fig. [2]). To improve on 
these deficiencies and to reach a higher confidence in the 
gaps of Refs. [ToLflll]. we generate a new parameterization 
of the Skyrme EDF, adding to the Lyon protocol [lij the 
constraint that both SNM and PNM HF effective masses 
shown in Fig.[5]are reproduced around saturation density. 
The refitted effective masses in both SNM and PNM are 
shown in Fig. [7l By default, the isoscalar effective mass 
reproduces well the microscopic one around saturation, 
whereas the isovector one is clearly improved compared 
to SLy4. However, the isoscalar effective mass is only 
marginally improved compared to SLy4 at low densities 
where it overestimates the HF result. In fact, a better 
reproduction of the low-density behavior of the effective 
mass in SNM requires extending the analytical form of 
the Skyrme functional [35j . This is underway, in partic- 
ular through the design of non-empirical Skyrme EDFs 
based on the development of the density-matrix expan- 
sion for low- momentum interactions [45l |46| . 

Pairing gaps calculated in SNM and PNM with the 
refitted Skyrme EDF and the soft NN interaction as a 
pairing kernel are presented in Fig. [8j With the improved 
isovector effective mass the gaps are more satisfactory in 
PNM compared to using SLy4 (see upper panel of Fig. Q|, 
but they are essentially unchanged in SNM, where in par- 
ticular, microscopic pairing gaps are overestimated for 
k F < 1.2 fin -1 . 



2. Pairing gaps in semi-magic nuclei 

Employing the refitted Skyrme EDF and the soft NN 
interaction as the pairing kernel, we compute neutron 
and proton pairing gaps in semi-magic nuclei by solving 
HFB equations in spherical symmetry Theoretical 
gaps are provided by Alcs which denotes the diagonal 
pairing matrix element Aj corresponding to the canoni- 



ing spectrum generated from Skyrme EDFs characterized by 
m* /m = 0.7 at nuclear saturation density. However, 40 Ca con- 
stitutes an anomaly, because the spectrum generated by such 
EDFs is unnaturally dense around the Fermi energy l42t U'.'ti . In 
addition, a meaningful microscopic calculation of single-particle 
spectra must include 3N interactions. 



cal single-particle state </>j whose quasi-particle energy 5 
is the lowest. Our results are presented in Fig. [5J Ex- 
perimental gaps extracted from binding energies through 
three-point mass differences centered on odd-mass nu- 
clei [47] are shown as a reference. The pairing gaps are 
essentially identical to those obtained using SLy4 [ll[ 
where data exist. The improved isovector effective mass 
leads to a tiny increase (decrease) of neutron (proton) 
gaps in neutron-rich nuclei. Note that, although the dis- 
crepancy between the refitted Skyrme effective mass and 
the microscopic results at low densities in SNM implies 
an uncertainty for the present results, we expect this to 
be small, as it is unlikely that such densities weigh sig- 
nificantly in pairing gaps of (non-halo) nuclei. 

Keeping in mind the necessity to confirm the effective- 
mass approximation through a systematic comparison of 
HF and Skyrme single-particle spectra in doubly-magic 
nuclei, one can conclude that the NN-only results of 
Refs. [13, [H| are presently put on a rather solid basis. 
The most important of these conclusions is that neutron 
and proton pairing gaps in semi-magic spherical nuclei 
are approximately accounted for using the partial 
wave of soft interactions at first order. This result is 
valid [ll[ over the broad cutoff range of A s» 1.8—3.0 fm _1 
that characterizes perturbative NN interactions 
This finding is somewhat puzzling as it indicates that 
neglected contributions, such as higher partial waves, 3N 
interactions and coupling to density, spin and isospin col- 
lective fluctuations for both the normal self-energy and 
the pairing interaction, may have a small net effect on 
pairing gaps in nuclei 6 . Such a conjecture needs to be 
confirmed by incorporating explicitly all neglected con- 
tributions in a consistent way. 



B. Hard interaction 

Let us now turn to the hard NN interaction case. The 
findings discussed in Section IIII B 21 are useful to assess 
the validity of the results of Refs. [HI, [l3j]. Indeed, it 
is natural to expect that the uncertainties in pairing 
gaps generated in INM by the momentum-independent 
effective-mass approximation propagate to finite nuclei. 
This is the case because the uncertainty relates to high- 
momentum modes that are of a similar nature in homoge- 
neous and non-homogenous systems. Practical applica- 
tions combining a hard NN interaction as pairing kernel 
with a momentum-independent effective mass provided 
by a phenomenological Skyrme functional would require a 
fine tuning of the momentum-independent effective mass 
and of the momentum-independent Z-factor when solv- 
ing HFB equations. Due to the energy dependence of 



J The acronym LCS stands for Lowest Canonical State. 
6 Of course, the neglected contributions do not have to be individ- 
ually small. 
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FIG. 9: (Color online) Neutron and proton LCS pairing gaps computed in semi-magic nuclei using a soft/hard NN interaction 
and the corresponding refitted Skyrme EDF whose effective masses are shown in Fig. [TJ/Fig. [TT] In the hard-interaction case, 
a Z-factor is taken into account when solving the HFB equations. Calculations include the Coulomb interaction in the proton 
pairing kernel [ill ]. 
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FIG. 10: (Color online) Neutron and proton LCS pairing gaps computed in semi-magic nuclei using a hard NN interaction and 
the corresponding refitted Skyrme EDF whose effective masses are shown in Fig. 1121 The ^-factor and e-mass are not taken 
into account when solving the HFB equations. Calculations include the Coulomb interaction in the proton pairing kernel |Tl| . 



the BHF normal self-energy, a consistent calculation of works, 
pairing gaps within the Brueckner expansion requires the 
inclusion of Z-factors. However, since the HFB calcula- 
tions of Refs. [HI, [TH are performed without Z-factor, 
both classes A^ ard and Aj arcl (see Section IHE)) are inves- 
tigated in this study in order to make contact with these 
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FIG. 11: (Color online) New Skyrme effective masses con- 
strained to reproduce BHF total effective masses obtained 
from the hard NN interaction. We have used a simple func- 
tional form to fit the Z-factors. 



1. Refitted Skyrme effective mass and Z -factor 

In order to assess the uncertainties of the pairing gaps 
in hnite nuclei due to the effective-mass approximation 
schemes, we use the INM results of Section IIII A 21 for 
the construction of new Skyrme functionals for hard in- 
teractions. We again consider the two classes of Table [Til 

A fard Wlth m * = m *pe/av and Z = Z P e/av as well aS A£ Md 

with m* = m£ and Z = 1. For the two classes. 

A, / Lb V ' 

the difference of the gaps using the point-evaluated and 
averaged quantities provides a range for the uncertainty 
in the pairing gaps. 



FIG. 12: (Color online) New Skyrme effective masses con- 
strained to reproduce BHF fc-masses obtained from the hard 
NN interaction. 



The refitted Skyrme effective masses in SNM and PNM 
are shown in Figs. QT] and [T2] As in the soft-interaction 
case, the effective mass is satisfactory around saturation 
density in both SNM and PNM. Only the fit to the point- 
evaluated fc-mass in SNM is somewhat problematic due 
to the density dependence of the Skyrme parametriza- 
tion. Reproducing well the BHF effective mass at low 
densities would require to extend the analytical form of 
the Skyrme functional (3f|. 



2. Pairing gaps in semi-magic nuclei 

Employing the refitted Skyrme EDFs and the hard NN 
interaction as the pairing kernel, we compute neutron and 
proton pairing gaps in semi-magic nuclei. The results are 
shown in Figs. [H andfTUl 

We observe a reduction of the pairing gaps obtained 
within the Aff ard class compared to those obtained from 
soft NN interactions. This is consistent with the results 
obtained in INM and discussed in Section [TTTB] This dif- 
ference is genuine and reflects that many-body expansion 
schemes depend on the resolution scale and that finite- 
order results for soft and hard interactions are in general 
not equivalent and immediately comparable. Effectively, 
the Z-factor accounting for the energy dependence of the 
BHF self-energy is largely responsible for the smallncss of 
the pairing gaps computed at lowest order in the hard NN 
interaction. Omitting the Z-factor and e-mass makes the 
present calculation formally similar to those performed 
in Refs. [12, EH which (using the SLy4 parameterization 
with an isoscalar effective mass of m* jm = 0.7 at satu- 
ration density) led to pairing gaps smaller by a factor of 
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two-thirds compared to those obtained with soft interac- 
tions [Tlj . Using the new Skyrme parametrization, such 
a calculation leads to larger gaps than in class Ajf ard , but 
still generally smaller than with soft interactions. How- 
ever, for both classes Aff ard and A^ ard , the dependence 
of the pairing gaps on the effective-mass approximation 
scheme is significant. For class A^ ard the uncertainty is 
w 0.5 MeV and for class A^ ard it is w 0.25 MeV, which 
both constitute a substantial fraction of the gap strength. 
This demonstrates that the effective-mass approximation 
problem is present irrespective of the inclusion or omis- 
sion of the energy dependence of the self-energy. 

V. CONCLUSIONS 

The present paper complements recent works directed 
towards the construction of non- emp irical energy func- 
tional for nuclei 0, HI E3, El Ef. We have studied 
neutron 1 So pairing gaps with special attention to the 
consistency of the pairing interaction and normal self- 
energy contributions. In nuclear matter, we calculated 
the normal and anomalous parts of the self-energy consis- 
tently at first order in the expansion scheme for soft and 
hard NN interactions. Our results also provide new con- 
straints to empirical Skyrme functionals. We have found 
that T = 1 pairing gaps obtained from low-momentum 
interactions depend only weakly on approximations to 
the normal self-energy, while gaps from hard potentials 



are very sensitive to the effective-mass approximation 
scheme. This is because a momentum-independent ef- 
fective mass does not approach the free mass at high 
momenta, but for hard interactions the high-momentum 
modes are not decoupled. The same conclusion has 
been reached for calculations of pairing gaps in finite nu- 
clei. This is problematic for hard NN interactions when 
employed in conjunction with standard empirical EDFs 
which are of low-momentum character. Although a com- 
plete first order calculation is needed, where the Skyrme 
EDF is replaced by a microscopic HF calculation includ- 
ing 3N forces, our results put the effective- mass approx- 
imation used in Refs. [T5L IlO ] on a rather solid basis. 
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